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A many-body theory approach is developed for the problem of positron-atom scattering and 
annihilation. Strong electron-positron correlations are included non-perturbatively through the 
calculation of the electron-positron vertex function. It corresponds to the sum of an infinite series 
of ladder diagrams, and describes the physical effect of virtual positronium formation. The vertex 
function is used to calculate the positron-atom correlation potential and nonlocal corrections to the 
electron-positron annihilation vertex. Numerically, we make use of B-spline basis sets, which ensures 
rapid convergence of the sums over intermediate states. We have also devised an extrapolation 
procedure that allows one to achieve convergence with respect to the number of intermediate-state 
orbital angular momenta included in the calculations. As a test, the present formalism is applied to 
positron scattering and annihilation on hydrogen, where it is exact. Our results agree with those of 
accurate variational calculations. We also examine in detail the properties of the large correlation 
corrections to the annihilation vertex. 
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' I. INTRODUCTION 

. ^ i The interaction of a low-energy positron with a many-electron atom is characterised by strong correlation effects. 
^ ■ Apart from the dynamic polarisation of the electron cloud by the field of the positron, the positron can also form 

positronium (Ps), by picking up one of the atomic electrons. When the positron energy is below the Ps-formation 
'"q^. threshold, eps = / + ii^is(Ps) = I — 6.8 eV, where / is the atomic ionisation potential, positronium formation is a 
I— —I' virtual process. Nevertheless, its role in the positron-atom interaction is very important The main aim of 

this work is to develop a many-body theory approach which accounts accurately for both correlation effects, and to 
" ' ' test it for positron scattering and annihilation on hydrogen. 

^ [ The study of positron interaction with matter is a topic of fundamental interest Q. Positrons have also found 
I ■ many useful applications. They are a very sensitive probe of the presence of defects in materials @. The recent 
[ development of a scanning positron microscope may lead to an increased use of positrons for quality control 
CO i: of materials, particularly in the semiconductor industry. In medicine, positron emission tomography, or PET, has 
■ become a standard means of medical imaging (see, e.g., 0). A proper understanding of how positrons interact with 
, matter at the fundamental level of atoms and molecules, is therefore essential. 

The interaction of low-energy positrons with atoms has presented a challenge to the theorist for many decades. 
^ , The study of positron scattering from atoms was initially seen as a useful complement to work on electron scattering, 
^ ■ particularly in helping to understand the role of the exchange interaction. However, although the exchange interaction 
^ ' is absent, it was quickly realised that the positron-atom problem is more complex than the electron case. The attractive 
induced polarisation potential tends to cancel and even overcome the static repulsion of the positron by the atom 
at low energies. The positron may also temporarily capture one of the atomic electrons in a process known as 
virtual positronium formation. The need to account for these effects requires an elaborate and accurate theoretical 
description. 

k> ', For small systems, such as hydrogen and helium, accurate results were obtained through the use of variational 
j_j ■ methods fiol ITlL [T3 . H^ . Positron scattering from alkali atoms which have a single valence electron, has 
I been calculated extensively using a coupled-channel method with pseudostates More recently, positron and Ps 
interaction with atoms with few active (valence) electrons has been studied using the stochastic variational method 
(SVM) and configuration-interaction-type approaches [Til ITtI ITM ITqI |20J . However, it is difhcult to extend these 
methods to larger atomic systems with many valence electrons, e.g., the noble gases. 

An attractive alternative to few-body methods is many-body theory [23|. It lends itself naturally to the study 
of problems where an extra particle interacts with a closed-shell target. The use of diagrams makes this method 
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both descriptive and intuitive, and allows one to take many-particle correlations into account in a systematic way. 
Many-body theory has been successful in the study of photoionisation ,2% . and in problems involving electrons, such 
as electron-atom scattering [23, 0, HE 11^ > negative ions |23, HE 1291 ISf l . and precise calculations of energies and 
transition amplitudes in heavy atoms with a single valence electron |3ll l32| . The application of many-body theory to 
low-energy positron interactions with atoms has met with more difficulty. 

Many-body theory utilises techniques originating in quantum field theory. It describes the terms of the perturbation 
series in the interaction between particles diagrammatically. The difficulty in applying this approach to the interaction 
of positrons with atoms arises from the need to take into account (virtual) Ps formation. Being a bound state, Ps 
cannot be described by a finite number of perturbation-theory terms. Hence, an infinite sequence of the 'ladder' 
diagrams must be summed. 

The first attempt to apply many-body theory to the positron-atom problem was by Amusia et al in 1976 0], 
who used a crude approximate method of accounting for virtual Ps formation in He. A better approximation for 
the virtual Ps- formation contribution was devised in [^, and applications to various atomic targets, including noble 
gases, were reported I^EISH^- In particular, a reasonable description was obtained for positron scattering from 
noble-gas atoms, which highlighted the presence of positron-atom virtual levels in Ar, Kr and Xe. On the other hand, 
application of the same approximation to positron- atom annihilation showed that it was clearly deficient. 

In spite of the approximate treatment of virtual Ps formation, many-body theory calculations for Mg, Cd, Zn 
and Hg were the first to provide credible evidence that positrons can bind to neutral atoms [s^. Two years later 
positron-atom binding was proved in a stochastic variational calculation for Li .36j . At present the list of atoms 
capable of binding positrons has expanded dramatically, SVM and configuration-interaction calculations confirming 
positron binding to Mg, Cd and Zn |33,|3E|- This topic is now of major interest in positron physics (see [ssl Isil l for 
useful reviews). 

In this paper, new techniques will be outlined that allow the exact calculation of the electron-positron ladder 
diagram sequence which accounts for virtual Ps formation. This approach enables many-body theory to provide 
accurate information on the elastic scattering, annihilation and binding of positrons on atoms and negative ions at 
energies below the Ps formation threshold. 



II. MANY-BODY THEORY METHOD 



A. Dyson equation 

A conventional treatment of positron scattering from an iV-electron target would start from the Schrodinger equation 
for the total wavefunction for the iV + 1 particles. In many-body theory we start instead from the Dyson equation 
(see, e.g., ^M). 

{Ho + I]e)V'e = £^e, (1) 

where ip^ is the single- particle (quasi-particle) wavefunction of the positron, s is its energy, and Hq is a central-field 
Hamiltonian of the zeroth approximation, which describes the motion of the positron in the static field of the target. 
The many-body dynamics in Eq. is represented by E^, a nonlocal energy-dependent correlation potential. This 
quantity, also known as the optical potential, is equal to the self-energy part of the single-particle Green's function of 
the positron in the presence of the atom |4lj | . Due to its nonlocal nature operates on the quasi-particle wavefunction 
as an integral operator, 

= J Ee(r,r')^e(r')dr'. (2) 

For hydrogen Hq may simply be taken as the Hamiltonian of the positron moving in the electrostatic field of the 
ground-state atom, Hq — — ^V^ -I- U{r), where U{r) = (1 -|- r~^)e~^^ ^3 ^^e atomic units throughout). For 
systems containing more than one electron the Hartree-Fock (HF) Hamiltonian (without exchange, for the positron) 
is the best choice. The correlation potential is given by an infinite perturbation series in powers of the residual 
electron-electron and electron-positron interaction. Inclusion of the electrostatic interaction in Hq and the use of 
the HF approximation for the target electrons means that the perturbation-theory expansion for starts with the 
2nd-order diagrams, and that the diagrams do not contain elements which describe the electrostatic potential l43j| . 

Owing to the spherical symmetry of the problem, Eq. can be solved separately for each positron partial wave. 
So, in practice one deals with radial quasiparticle wavefunctions, P^i{r), related to ipe by V'e(r) — i"^^Pei{''^)Yim{fl), 
where YimiS^) is the spherical harmonic for the orbital angular momentum I. Accordingly, the self-energy operator is 
also found for each partial wave separately, as Ee''(r, r'), see Eq. (|A1|I in the Appendix. 
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B. Correlation potential 



Figure ^ shows the lowest-order terms of the diagrammatic expansion for the correlation potential E, or more 
precisely, for the matrix element (e'jE^le) of the correlation potential calculated at some energy E between the 
positron states e and e' . The leading 2nd-order diagram, (a) in Fig. ^ corresponds to the following expression. 



{e'n\V\fii^){i^^\V\ne) 
E - - Ea + £n + iO' 



(3) 



where V is the electron-positron Coulomb interaction, the sum runs over all intermediate positron states v, excited 
electron states fj, and hole states n, and iO is an infinitesimal positive imaginary quantity. 




+ ... 



FIG. 1: Diagrammatic expansion of the positron-atom correlation potential E. Shown are the 2nd-order and main 3rd-order 
contributions. The top line in all the diagrams corresponds to the positron. Other lines with the arrows to the right describe 
excited electron states, while those with the arrows to the left correspond to holes, i.e., the electron states occupied in the atomic 
ground state. Wavy lines are the electron-positron or electron-electron Coulomb interactions. 



It is easy to see from Eq. (O that at low energies E the diagonal matrix element (e|E^ |e) is negative. This means 
that the 2nd-order contribution to the correlation potential. Fig. Q^a), describes attraction. In fact, this diagram 
accounts for the main correlation effect in low-energy scattering, namely the polarisation of the atom by the charged 
projectile. At large distances it leads to a well-known local polarisation potential, 

EgHr,r')^-^<5(r-r'), (4) 

where a is the static dipole polarisability of the atom in the HF approximation, 

Besides the 2nd-order diagram. Fig. QJshows the main 3rd-order contributions. A complete list of 3rd-order diagrams 
includes mirror images of the diagrams (f) and (g). There are also a few more diagrams similar to diagram (h), where 
the positron line is connected to the atomic excitation part by a single line. They represent correlation corrections to 
the HF electron charge density of the ground-state atom. Such corrections are much smaller than other correlations 
effects jS^, and can be neglected. The total number of the 3rd-order diagrams in the positron-atom problem is 
considerably smaller than that in the electron case (see, e.g., [23 )i where one needs to allow for the exchange between 
the incident and core electrons. 

Comparing diagrams (c), (e), (f) and (g) with (a) in Fig. ^ we see that they represent corrections to the leading 
polarisation diagram, S^^^ due to electron correlations within the atom. Interaction between electron-hole pairs can 
in principle be included in all orders, which would correspond to the random phase approximation (RPA) treatment 
of atomic polarisation filj. On the other hand, if the two hole orbitals in diagrams (c) and (e) are identical, these 
diagrams together with similar higher-order terms, are easily incorporated within the 2nd-order diagram by calculating 
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the excited electron states fi in the field of the atom with a hole in this orbital [24l |. These approximations, and even 
the "bare" 2nd-order approximation (with exchange diag rams added in both cases), give good results in electron-atom 
scattering and negative ion problems IIJ, I29t l8(il| . 

However, for the positron-atom problem the approximation based on diagrams (a) and corrections of types (c), (e), 
(f), and (g), proved to be deficient ^^tlj. In addition one must include the diagram Fig. H (b) and higher-order 
diagrams in which the positron interacts with the excited electron in the intermediate state. Fig. |21 This sequence 
of diagrams accounts for virtual Ps formation. It is important that it is summed to all orders, since in quantum 
mechanics a bound state (here, Ps) which is absent in the zeroth approximation cannot be described by a finite 
number of perturbation theory terms. 



e Vj V2 e' s Vj V2 e' 




FIG. 2: Virtual Ps contribution to the positron-atom correlation potential S. 



Summation of the diagrammatic sequence shown in Fig. |2] is done by calculating the electron-positron vertex 
function P, defined here as the sum of the electron-positron ladder diagrams. Fig. |3 and denoted in the diagram by 
the shaded block. 



Vi V2 Vi V2 Vi V V2 




1^1 1^2 ^1 1^2 1^1 P-2 



FIG. 3: Electron-positron ladder diagram sequence and its sum, the vertex function F. 



Comparing the left- and right-hand sides of the diagrammatic equation in Fig. |31 we see that the vertex function 
satisfies the following linear equation, 

(i/2M2|Pe|M1'^i) = (i^2M2|I^|Mlf^l) + 2^ —-7: ■ (6) 

The vertex function depends on the energy E. It has the meaning of the electron-positron scattering amplitude in 
the field of the atom. In the lowest-order approximation Te — V. 

Once the vertex function is found, the contribution of virtual Ps to the correlation potential. Fig. |21 is obtained as 

Since the electron-positron Coulomb interaction is attractive, ^ < 0, the terms of the perturbation series in Fig. El 
have the same sign. This explains why their sum, and hence, the contribution of virtual Ps formation to the positron- 
atom attraction, are large. Note that for electron scattering {V > 0) this series is alternating. As a result the net 
contribution of the diagrammatic series on the left-hand side in Fig. [3 is small, and its omission in the electron- atom 
correlation potential does not give rise to large errors. 

The Ps-formation contribution to the correlation potential was previously approximated by using the free Is-state 
Ps propagator orthogonalised to the ground-state electron wavefunctions [j, yj] , 

/./|..(^|,^ ^, f {£'n\V\i',s^K){^is,K\V\ne) dK 

where xE'is.k = {Sn)-'^/'^ ex.p{-\r - r'|/2)exp[iK • (r -I- r')/2] is the wave function of Ps(ls) with momentum K, 
Els + K'^/4: is the energy of this state, and the tilde above ^'is,k in Eq- © indicates orthogonalisation. This 
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approximation is suitable for positron scattering from the targets where ground-state Ps formation dominates, e.g., 
hydrogen or noble gas atoms. It also allows one to consider positron scattering above the Ps-formation threshold, 
where the correlation potential acquires an imaginary part due to the pole in the integral in Eq. © QB| ■ At the 
same time, the ground-state Ps propagator fails to describe the short-range electron-positron correlations crucial for 
the calculation of the annihilation rates "34'| . By the uncertainty principle, small separations invoke contributions of 
highly excited states of the Ps internal motion, not included in the Ps(ls) propagator. In contrast, our present method 
based on the summation of the ladder diagram series, Eqs. © and {71), is consistent and complete. It accounts for 
all (virtual) intermediate states of the electron-positron pair. 

For positron scattering on hydrogen, only a few types of diagrams contribute to S^;, since only one hole can be 
created. Moreover, the interaction of the intermediate-state electron and positron with the hole [diagrams (c) and 
(d) in Fig. ^ can be taken into account by calculating the intermediate electron and positron wavefunctions in the 
Coulomb field of the nucleus. In this case, the correlation potential is given by the sum of the 2nd-order diagram and 
the virtual Ps contribution. Figs. H^a) and|2 and 

Ss = 4'^+4'^' (9) 

is the exact correlation potential. In particular, the long-range asymptotic behaviour of S^'' at low energies, Eq. Q), 
contains the exact polarisability of hydrogen, a = |. 

C. Scattering 

Rather than solving the Dyson equation for the quasiparticle wavefunction in the coordinate representation, it is 
easier to work with the self-energy matrix, 

(e'|S£|£) = J ^;,{v)i:E{r,r')^e{r')drdr', (10) 

where ipi^ are the positron eigenfunctions of the HE (or ground-state hydrogen) Hamiltonian Hq, 

HQipe=£Ve, (11) 

with a given angular momentum /, (y5e(r) = P^i{r)Yimi^)- Since the static potential of the atom is repulsive, all 
positron states (p^ lie in the continuum (e > 0). The radial wavefunctions are normalised to a J- function of energy in 
Rydberg, S{k^ — k' ), where k is the positron momentum. This corresponds to the asymptotic behaviour 



i{r) ~ (7rfc)-i/^ sin (fcr - ?7r/2 + Sl°^^ , (12) 



where S^^^ is the scattering phase shift in the static potential. 

The matrix can be used to obtain the phaseshifts directly [2^. First, a 'reducible' self-energy matrix (e'lEsle) 
is found via the integral equation, 

{e \^e\£) = {£ \^e\£) +7^ E - e" ^ ^ 

where V means the principal value of the integral. The phase shift is then given by, 

Si = Sl''^ +ASi, (14) 

where 

tanA(5; = -27r(e|I]e|e), (15) 

determines the additional phaseshift ASi{k) due to correlations, at the positron energy e. 

Once the reducible self-energy matrix has been found, the quasiparticle wavefunction (also known as the Dyson 
orbital) can be found via, 

V^^r) ^iPeir)+V f ^,,(r)iffl£^d£'. (16) 
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In order to normalise the quasiparticle radial wavefunction at large distances to 

P.iir) - (7rfc)-i/2 sin (kr - Itt/2 + S^"^ + ASi^ , (17) 
the function obtained from the right-hand side of Eq. I|16f) must be multiplied by the factor, 



cos ASi 



1+ 2^(e|Ee|e) 



-1/2 



(18) 



D. Positron annihilation 

The annihilation rate A, of a positron in a gas of number density n is usually expressed in terms of the effective 
number of electrons, Z^b, which contribute to annihilation on an atom or molecule |4fil |. 

A = irr'^mZcS, (19) 

where ro is the classical electron radius and c is the speed of light. Equation (|19|l defines Zcs as the ratio of the 
positron two-photon annihilation cross section of the atom to the spin-averaged two-photon annihilation cross section 
of a free electron in the Born approximation [4^ . Annihilation takes place at very small electron-positron separations, 
h/{mc) ~ 10~^ a.u. Hence, for nonrelativistic positrons it is determined by the electron density at the positron, and 
Zeff can be calculated as |4m |. 



N 
i=l 



j |*(ri, r2, . . . , TAT, r)|2 5{v, - r)dri . . . drwdr, (20) 



where ^'(ri, r2, . . . , r^r, r) is the full [N + l)-particle wavefunction of the N electron coordinates and positron 
coordinate r. The wavefunction is normalised to a positron plane wave at large positron-atom separations, 

4'(ri,r2,...,rAr,r) ~$o(rl,r2,...,r^r)e*•^ (21) 

where $o(ri, r2, ■ ■ ■ , '^n) is the atomic ground-state wavefunction, and k is the incident positron momentum. 

Although Zcff is basically a cross section, Eq. (l2Uf) has the form of a transition amplitude. This enables one to apply 
the apparatus of many-body theory to this quantity directly 0, |3 . In this "transition amplitude" the initial and 
final states are identical, and the electron-positron two-body operator, 5{Yi — r), plays the role of a perturbation. 
The positron energy in the initial and final states is the same, e — A:^/2, and (owing to the spherical symmetry 
of the target) the perturbation conserves the positron angular momentum I. Therefore, the contribution of each 
positron partial wave to ZcB can be determined separately. The corresponding many-body diagrammatic expansion 
is presented in Fig. 0] 

The analytical expression for the Oth-order diagram. Fig. Efa), is 

^cff = E / V^:«^:(ri)<5(r - ri)^„(ri)V'e(r)drdri = ^ / |^„(r)n^,(r)pdr, (22) 

n n 

where i/'e is the positron wavefunction, tpn is the wavefunction of the hole, and the sum over n runs over all holes, 
i.e. orbitals occupied in the target ground state. This contribution is simply an overlap of the electron and positron 
densities, ^„ |(/9„(r)p and |i/'e(r)p, respectively. 

The two Ist-order "corrections", Fig. 01b) and (c), are identical, and their contribution is 



= 2 E ^'^'^l l^^^^''^! I'^^^ (23) 

cf. Eq. |(2Jl for the 2nd-order contribution to the correlation potential. In the calculations of annihilation we assume 
that the positron energy is below all other inelastic thresholds, hence, we have dropped «0 in the energy denominator. 
Physically, the Ist-order diagram describes positron annihilation with an electron "pulled out" from the atom by 
the positron's Coulomb field. Calculations in Refs. 0, 13 showed that for noble-gas atoms the size of the Ist-order 
corrections is approximately equal to the Oth-order contribution. This means that higher-order terms must also be 
taken into account. 
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(a) (b) (c) (d) 




(e) (f) (g) (h) 

FIG. 4: Diagrammatic expansion of Zctt. Apart from the 0th and 1st order diagrams (a), (b) and (c), the main types of 
2nd-order diagrams are shown. The external lines of these diagrams (e) represent the wavefunction of the incident positron. 
The solid circle corresponds to the electron-positron 5-function annihilation vertex. 

Diagrams (d)~(h) in Fig. ^illustrate the main types of 2nd-order corrections to the annihilation vertex. It is also 
important to consider corrections to the incident positron wavefunctions denoted by e. However, these corrections are 
included in all orders in the positron quasiparticle wavefunction, obtained from the Dyson equation , or via equation 
(|16|l . Hence their contribution to Zcs is accounted for by using the positron Dyson orbitals ip^ in the calculation of 
the annihilation diagrams. 

The presence of the (5-function operator in the annihilation diagrams enhances the importance of small electron- 
positron separations. For this reason, the most important diagrams in ZeS are those with the Coulomb interactions 
between the annihilating pair, e.g., the 2nd-order diagrams (d) and (e) in Fig. ^ (the latter together with its mirror 
image), and similar higher-order terms. A complete all-order calculation of their contribution is achieved by using 
the vertex function, as shown in Fig. [S] Note that for hydrogen this set of diagrams is exhaustive, provided the 
intermediate electron and positron states are calculated in the field of the bare nucleus. Previously diagrams (c)~(f) 
in Fig. 13 have only been estimated We will see that the ability to calculate these diagrams accurately is crucial 
for obtaining correct values of Zcb- We will also see in Sec. IIVI that the role of the vertex function (representing 
virtual Ps) in annihilation is much greater than in scattering. 
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(a) (b) (c) (d) 




(e) (f) 

FIG. 5: Main contributions to the positron-atom annihilation parameter Zcg. To account for the mirror images of the diagrams 
(b), (d) and (e), their contributions are multiplied by 2. 

Of course, for many-electron targets one can also consider other diagrams, e.g., (f)~-(h) in Fig. 0] In particular, 
diagram (f) describes screening of the electron-positron Coulomb interaction by other electrons [cf. Fig. (e)]. 
Diagram (g) can be viewed as the lowest-order "pick-off" annihilation contribution. Here the positron excites an 
electron-hole pair (a precursor of virtual Ps formation) and annihilates with an electron from one of the ground-state 
orbitals. Diagram (h) is independent of the positron energy. It represents one of the corrections to the HF ground-state 
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electron density, cf. Fig. ^h). Unlike the diagrams in Fig. [SJ these contributions are not systematically enhanced by 
the electron-positron Coulomb interaction at small distances. 

It is clear from Fig. 0]that most correlation corrections to the annihilation vertex, including the dominant sequence 
of diagrams in Fig. |S1 are nonlocal. As a result, the total Zcs can be written as 



Zes = Jj2\^nir)\'mr)\^dr + J V;(r)A,(r, r')Vs(r')drdr', 



(24) 



where Ae(r, r') represents the nonlocal correlation correction to the annihilation vertex. In the approximation of Fig. 
|31 it is equal to the sum of all diagrams (b)-(f) with the external positron hues e detached. Diagram (h) in Fig. ^ 
and similar diagrams which represent corrections to the HF electron density, could be included by replacing of the 
HF electron density |(^„(r)p in Eq. (|24|l with the exact target electron density Pe(r). Given the high accuracy of 
the HF density, this would make only a small change in Z^s- The structure of Eq. (|24|l shows that even when one 
uses the best single-particle positron wavefunction, the annihilation rate is not reduced to a simple (local) overlap of 
the electron and positron densities. 

There is an important physical difference between the correlation effects in positron scattering and annihilation. 
The key role played by the long-range polarisation potential for low-energy positrons means that large distances are 
important. Polarisation also emphasises the role played by the dipolc part of the positron-target Coulomb interaction 
and dipole target excitations. The contribution of virtual Ps formation to S^; is a short-range effect. The typical 
distances here are comparable to the radius of the atom, or the radius of ground-state Ps. The net effect of the 
strong positron-atom attraction brings about low-lying virtual s-states (see, e.g., ^3) Kr and Xe [s^l, or 

positron-atom weakly bound states, e.g., in Mg |3a. ISSj. In both cases, the positron scattering phaseshifts and the 
positron quasiparticle wavefunction in the vicinity of the atom vary rapidly as functions of the positron energy. 

On the contrary, in the annihilation diagrams the (5-function vertex emphasises small electron-positron separations. 
By the uncertainty principle, such small separations correspond to high-energy excitations in the intermediate states 
in the diagrams (b)-(f). Fig. [S| As a result, the non-local correction to the annihilation vertex, A^, has a weak energy 
dependence, and the energy dependence of the 2nd term in Eq. H24f) is almost entirely due to that of the positron 
wavef unctions. The only exception is when the positron energy approaches the Ps-formation threshold from below. 
Here the virtual-Ps contribution. Fig. I^Jf), rises sharply. Details of the threshold behaviour of Zcs are discussed in 
Ref. J8]. 

The large difference in the energy scales characteristic of positron scattering and annihilation has another physically 
important consequence. It turns out that the relative size of the annihilation vertex corrections in Z^s, i.e. the ratio of 
the 2nd term in Eq. I|24() to Z^^ , is about the same, whether are the positron wavefunctions in the repulsive static 
potential, ip^ = y:>e, or the Dyson orbitals which fully account for the positron-atom correlation potential. Numerical 
illustrations of this effect will be provided in Sec. IIVI 

Finally, we should mention that the correct normalisation of Z^s for positrons with angular momentum I is obtained 
by multiplying the diagrams in Figs. 0]and|Slby the extra numerical factor, 

47r2 

— (2^ + 1). (25) 

This follows from the structure of the positron wavefunction ipi^, which has the asymptotic behaviour of a plane wave 
e^^ '^ at large distances [cf. Eq. l|2J], 

. I — oo m—l 

= -yJlT. E ^'e^''Pei{r)YC^{n^)Yira{nr). (26) 
' ;=o m=-l 

To derive Eq. H25|) . one can use "01: as the external positron lines in an annihilation diagram, and perform averaging 
over the directions of k. 



III. NUMERICAL IMPLEMENTATION 

A. Use of B-splines and convergence 

To evaluate the diagrams of the correlation potential E and annihilation parameter Z^b, one first needs to generate 
sets of electron and positron HF basis states. These are then used to calculate matrix elements of the Coulomb 
and 5-function operators, the main building blocks of the diagrams. Evaluation of the diagrams requires summation 
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over complete sets of electron and positron intermediate states, including integration over the electron and positron 
continua. 

To perform a numerical calculation, the continuous spectrum can be discretised. The simplest way of doing this is 
by placing the system in a spherical cavity of radius R. Setting the wavefunctions to zero at the boundary will result 
in a discrete spectrum of eigenstates with an approximately constant stepsize in momentum space, 

A/c « tt/R. (27) 

If the value of R is sufficiently large (i? 3> i?at, where i?at ^ 1 a.u. is the size of the atom), the presence of 
the boundary will not affect the quantities calculated. Indeed, for the positron energy below all inelastic thresholds 
(except, of course, annihilation), the intermediate states in the diagrams are virtual^ and no particle in an intermediate 
state can escape to infinity. 

The drawback of this procedure is that for a suitably large i?, the stepsize in momentum is small, e.g., for i? = 30 a.u., 
Afc « 0.1 a.u. Hence, one would need large numbers of intermediate states to achieve convergence. Note that the 
actual upper energy limit depends on the quantity in question. Thus, diagrams in Z^g converge more slowly than 
those of the correlation potential S, because of the greater role of small electron-positron separations and high orbital 
angular momenta in the former. However, as a rough guide, summing up to the energy of 10^ a.u. should be sufficient 
for both. The question of the number of intermediate states is especially important for the calculation of the vertex 
function F, which is a Nr x Nr matrix [see Eq. ©], where Nr ^ N'^Hmax + 1), being the number of electron or 
positron states in each partial wave and /max being the largest orbital angular momentum included. It is clear that 
here the simple cavity quantisation cannot work. 

Instead, to achieve an accurate and economical span of the continuum we use B-splines |49| . B-splines of order k are 
n piecewize polynomials of degree k — 1 defined by a knot sequence rj which divides the interval [0, R] into n — k + 1 
segments |5Qj. The basis states are obtained by expanding the radial wavefunctions Pi{r) in terms of B-splines Bi{r), 

Pi{r)=Y,cf''B,{r), (28) 

and finding the eigenvectors and eigenvalues of the radial part of the HF (or hydrogen atom) Hamiltonian for each 
orbital angular momentum I by solving the generalised eigenvalue problem, 

Y^H^.cf^ =eY,Q,,cf\ (29) 

where Hij = {Bi\H'^\Bj), and Qij — {Bj^\Bj). Prior to solving Eq. (|29|l . the ground-state atom HF Hamiltonian is 
generated by a conventional HF routine (5lj. Note that in the sums over i and j in Eqs. (|28|) and H29|l the first and 
last splines are discarded to implement the boundary condition Pi(0) = Pi{R) = 0, leaving one with a set of n — 2 
eigenstates for each electron and positron orbital angular momentum. When Eq. H29|l is solved for the electron, it 
yields the wavefunctions of the orbitals occupied in the atomic ground state (holes), as well as those of the excited 
states (particles) . The exact energies of the excited electron and positron states are determined by the B-spline radial 
knot sequence. 

It is instructive to try to design an ideal distribution of energies of a discrete set spanning the continuum. Qual- 
itatively, at low energies (e <C 1 a.u.) the continuous spectrum states oscillate slowly, and the contribution of large 
distances in the matrix elements is important (hence, the need for a large R). As the energy of the states increases, 
the range of important distances becomes smaller and smaller. Indeed, the matrix elements then contain rapidly 
oscillating factors of e**''' type, which means that the dominant contribution comes from r ^ k~^. Therefore, one does 
not need a large value of the cavity radius R for the higher-energy states. More specifically, one can estimate the 
necessary radius as i? ~ a/k, where a is a number greater than unity. Combining this with the cavity quantisation 
condition (|27|) . one obtains Ak/k ~ 7r/a, which yields the following grid in momentum space: 

fcj. = fcoe''^ (30) 

where kg is the lowest momentum, /3 = 7r/a, and j = 0, 1, 2, . . .. Thus, it appears that the optimal momentum and 
energy grids are exponential. By choosing a small initial momentum fco <C 1 a.u. and /3 < 1 one ensures that the 
stepsize in momentum, Afc w is sufficiently small, to describe accurately the energy variation of the quantities 
summed. 

It turns out that basis sets generated by Eqs. (|28|l and (|29|) using an exponential radial knot sequence are a 
very close realisation of the exponential energy grid, and are effectively complete. In the present work we use n — AQ 
B-splines of order /c = 6 with a knot sequence 



=p(e-^~l). 



(31) 
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where p = 10~^ a.u., and a is determined by the condition rn-k+i = R- Figure shows the positive eigenenergies of 
the electron and positron basis states with I — for the hydrogen atom. Their distribution does indeed correspond to 
the exponential ansatz, Eq. with (3 k, a. The highest energies in the sets are about 10* a.u. This value is close 

to the magnitude of (pcr)"^, since the knot point closest to the origin determines the most rapidly varying eigenstate 
(by the uncertainty principle). In fact, it may not be necessary to include all 38 basis states in each partial wave in 
the many-body theory sums. In the calculations reported in this paper we use only about 15 lowest states, which 
span the energy range from threshold to 10^ a.u. 
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FIG. 6: Energies of the electron (open circles) and positron (solid circle) s-wave B-spline basis states in the field of the hydrogen 
nucleus, obtained using R — 30, k — 6 and n = 40. The first four electron states with negative energies (—0.50000, —0.12500, 
-0.05542 and -0.00246 a.u.) are not shown. 

Note that B-spline basis sets are used widely in atomic physics [i^, and that there are other basis sets which 
show a near-exponential spanning of the continuum. In particular, Laguerre basis states provide rapid convergence 
in close-coupling electron-atom scattering calculations (52 | . 



B. Calculation of the self-energy and annihilation diagrams 



The self-energy and annihilation diagrams are calculated by summation over the B-spline basis states, and the 
vertex function is found by matrix inversion from Eq. ©. The angular parts of the states are separated in the matrix 
elements and the angular variables are integrated over analytically. The actual expressions for the diagrams are given 
in the Appendix. The self-energy matrix and the vertex function are energy dependent. In practice, the self-energy 
matrix has been calculated at 8 energies spaced evenly from zero to the Ps formation threshold. Interpolation onto 
any required energy E is then used. 

Apart from the B-spline basis states, we also consider true positron continuum HE states (|12() . They are needed 
to evaluate the matrix elements (ejE^jje') and obtain the phaseshifts via Eqs. (|13|I - H15|) . Here we use 201 states that 
form an equidistant mesh in positron momenta of size AA: = 0.02. A transformation of the B-spline basis matrix 
elements (z|E£;|j) into (e|S£;|e') could be done using the effective completeness of the B-spline states on the interval 
[0,i?], 

{s\J:E\e')^Y.'^e\^){^\J:E\J){J\e'), (32) 

where {e\i) is the overlap of the HE state with the B-spline basis state. However, unlike the B-spline states which 
satisfy the zero boundary condition at r = _R, the continuous spectrum state P^i is finite at the boundary. To fix this 
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problem we insert a radial weighting function /(r) = R — r into equation H32() . which now reads 

(elE^ls'} = ^(e|/|z)(z|/-iE^rilj)0-|/|£'), (33) 

and calculate the "weighted" self-energy matrix (i|/~^E£;/~^|j), rather than (ijEslj). The singularity of /""'^ at 
r ~ R does not cause a problem, since the B-spline basis states in the Coulomb matrix elements involved, vanish at 
r ^ R. The same trick is applied in the calculation of the annihilation diagrams. 

To calculate (e|E£;|£') more accurately at low positron energies, where distances beyond r = R can be important, 
we make use of the long-range asymptotic form of the correlation potential l@J. The contribution oi r > R can be 
evaluated as 

/>oo 

/ Pei{r)[-^)Pe'i{r)dr, (34) 
with the correct value of the dipole polarisability a, and added to (elE^jje'}. 



C. Convergence with respect to the orbital angular momenta 

The use of a B-spline basis means that fast convergence is achieved with respect to the number of states with a 
particular orbital angular momentum. However, this leaves open the question of convergence with respect to the 
maximal orbital angular momentum of the electron and positron intermediate states included in the calculation. It 
has been known for a while that calculations of positron-atom scattering converge slowly with respect to the number 
of target aiigular momenta included in the expansion of the total wavefunction, notably slower than in the electron- 
atom case |53j. This is also true for the configuration-interaction-type calculations of positron-atom bound states 
[TtL IT^ Is^ and scattering |23| • Calculations of annihilation rates converge even more slowly 0, 0, . Physically, 
the slow convergence rate arises from the need to describe virtual Ps localised outside the atom by an expansion in 
terms of single-particle orbitals centred on the nucleus. 

The problem of the convergence rate with respect to the maximal orbital angular momentum has been investigated 
by the authors in Ref . ^55*1 . Using a perturbation-theory approach and the original ideas of Schwartz "56] , we derived 
asymptotic formulae that describe the convergence of the scattering amplitudes, or the phase shifts, and annihilation 
rates, or Zcs- The contribution of high orbital angular momenta probes small particle separations in the system. The 
difference between the convergence rates of the scattering and annihilation parameters is due to the presence of either 
the Coulomb interaction or the (5-function annihilation operator in the relevant amplitudes. 

The increments to the phaseshifts and Z^s upon increasing the maximum orbital angular momentum from A — 1 
to A, were found to behave as (A -|- 1/2)"'* and (A -I- 1/2) respectively. This means that if a series of calculations 
is stopped at some maximal angular momentum A = /max, the values obtained approach the ultimate (/max oo) 
values as follows: 



B 

(^max"Tl/2) ■ 



Z]^--^ ^ Zeff - (36) 



where A and B are some constants. They are determined together with 5i and Zcn by fitting Eqs. (|35|l and H36|) 
to the numerical data obtained for a range of /max- This extrapolation to /max — > oo is performed at each positron 
momentum value considered, and is especially important for obtaining correct values of Z^^. 



IV. RESULTS: SCATTERING AND ANNIHILATION ON HYDROGEN 



The theory outlined above can be readily applied to any closed-shell atom or ion. In this paper we would like 
to test it for the simplest possible target, the hydrogen atom. Since it contains only one electron, the correlation 
potential from Eq. (j^ and the Zcs diagrams in Fig. [Slgive an exact solution of the elastic scattering and annihilation 
problems, provided the intermediate electron and positron states are calculated in the field of H+. The key advance of 
the present many-body theory of positron-atom interactions relates to the calculation of the electron-positron vertex 
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function T, Eq. © and its incorporation in the self-energy and annihilation diagrams. It is mainly these features of 
the theory that a positron-hydrogen calculation is intended to test. 

In the numerical implementation we first generate the electron and positron B-spline basis sets. They are then used 
to evaluate the matrix elements, find the vertex function and calculate the self-energy and annihilation diagrams (see 
Appendix). Using the self-energy matrix, the phaseshifts are obtained by means of Eqs. (|13|I -H15 |) . and the positron 
Dyson orbital is calculated from Eqs. H16(l and (|18(l . In the end, the Dyson orbitals replace the positron HF states 
in the external lines of the annihilation diagrams, and final values of Zcft are obtained. To test the stability of the 
results with respect to the cavity radius, the calculations were performed with R — 15 and 30 a.u. 

To extrapolate the scattering phaseshifts and ZeS to l^ax ^ oo, as per Eqs. and the diagrams are 

evaluated for a range of maximal orbital angular momenta, Imax — 7-10. This procedure is illustrated by Fig. Qfor 
the phaseshifts and Fig. |Slfor Z^s, for the s-, p- and d-wave incident positron with momenta k = 0.2, 0.4, and 0.6 a.u. 
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FIG. 7: Convergence of the s-, p- and d-wave positron scattering phaseshifts on hydrogen with respect to the maximal orbital 
angular momentum Zmax for _R = 15 a.u. Circles, k — 0.2 a.u.; squares, k = 0.4 a.u.; diamonds, k = 0.6 a.u. 

Figures [3 and IHl show that the calculations have converged to the regime in which the asymptotic formulae H35|l for 
6i and (|36|) for Zcs, may be applied. The graphs also illustrate the point that the inclusion of high orbital angular 
momenta and extrapolation to Zmax — > oo is much more important in the calculations of annihilation, compared with 
scattering. For the positron momenta and partial waves shown, between 15% and 30% of the final value of Zcs is due 
to such extrapolation. Quantitatively, this contribution can be characterised by the ratio B/Z^s, see Eq. H36f) . given 
in Table Its increase with the positron angular momentum may be related to the greater role of the correlation 
corrections to the annihilation vertex in higher positron partial partial waves (see below). 

Figure 1^1 shows the s-, p- and d-wave phaseshifts for the total correlation potential @. They are in very good 
agreement with those from an accurate variational calculation (see Ref. IT3|), the discrepancy being of order 10"'^ rad. 
The values obtained with i? = 15 and i? = 30 a.u. are almost indistinguishable, except at low positron momenta. 
Here the results for i? = 30 are superior to those for R— 15. The larger cavity size allows for a better account of the 
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FIG. 8: Convergence of the s-, p- and d-wave contributions to Zcff for positron annihilation on hydrogen with respect to the 
maximal orbital angular momentum imax for ii = 15 a.u. Circles, k — 0.2 a.u.; squares, k — 0.4 a.u.; diamonds, k — 0.6 a.u. 

TABLE I: Values of B/Zcs which characterise the dependence of calculated Z^b on /max- 



Momentum 




Partial wave 




(a.u.) 


s 


P 


d 


0.2 


1.82 


2.63 


3.33 


0.4 


1.62 


2.58 


3.32 


0.6 


1.41 


2.51 


3.37 



long range —a/2r^ tail in the polarisation potential. 

Examining the phaseshifts allows us to compare the relative sizes of the polarisation and Ps-formation contributions 
to the correlation potential lO- The static positron-atom potential is repulsive, resulting in negative values of the 
phase shifts (thin solid curves). The inclusion of S, i.e., correlations, makes the low-energy phase shifts positive. 
Dashed curves in Fig. El show the phaseshifts obtained by including only the 2nd-order diagram E*^^^ (polarisation), 
while dot-dashed curves are those obtained with S''"' alone (virtual Ps-formation). We see that none of these results 
is close to the phaseshift obtained with the full S. This means that neither contribution dominates the correlation 
potential, and the inclusion of both polarisation and virtual Ps-formation effects is essential for solving the positron- 
atom problem. Of course, any calculation which produces accurate positron-hydrogen phaseshifts contains these 
contributions implicitly. The advantage of the many-body theory approach is that one can separate them, and get a 
better insight into the physics of the system. 
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FIG. 9: Positron-hydrogen s, p and d-wave scattering phaseshifts: long-dashed curve, many-body theory {R = 15 a.u.); sohd 
curve, many-body theory {R = 30 a.u.); circles, variational calculation 13]. Thin solid curve, static approximation; dashed 
curve, E'^'; dot-dashed curve, E'-'"'; dotted curve, E''"-' obtained with T = V (all R = 15). 

To illustrate the non-perturbative nature of virtual Ps formation we have also performed calculations that include 
the vertex function only to 1st order, Te = V, in S'-'") (dotted curves in Fig. inj. This approximation accounts only 
for about 50% of the total vertex function contribution. Note that the higher-order terms in T become even more 
important close to the Ps formation threshold {k « 0.7) in p and d waves. This is related to the virtual Ps becoming 
more "real" close to the threshold. 

We now turn to positron annihilation. Having solved the scattering problem accurately with the full S, we are now 
in possession of the best (quasiparticle) positron wavefunction, the Dyson orbital. Before using it in all annihilation 
diagrams, let us first look at the effect of the Dyson orbital on the Oth-order diagram. Fig. IS^a), for the s-wave 
Zgff. Figure [TUl shows that the Oth-order contribution, Eq. evaluated with the positron wavefunction in the 

static atomic potential gives values up to 20 times smaller than the accurate variational results [sJl ■ This situation is 
similar to that in positron annihilation on noble-gas atoms, where Eq. 1)22(1 evaluated with the static (HE) positron 
wavefunction underestimates experimental Z^s by a factor of 10^-10'^ at low positron energies 3, '34'|. It is natural 
that the use of the Dyson orbital, which is "aware" of the positron-atom attraction, in Eq. ((22(1 , leads to increased 
Zcs, and introduces a correct momentum dependence (dashed curve in Fig. I10() . This latter fact is in agreement with 
the general understanding of the origins of the energy dependence and enhancement of Z^s at low energies, and their 
relation to positron-atom virtual states ^, 34, 58] . However, the absolute values obtained are still about 5 times lower 
than the benchmark. 

The remaining 80% come from the nonlocal corrections to the annihilation vertex, diagrams (b)-(f) in Fig. |S1 The 
contributions of all the diagrams evaluated using the positron Dyson orbitals, and the total Z^g are shown in Fig. 1111 
for the positron s, p and d partial waves. 
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FIG. 10: Contribution of the 0th order diagram, Fig. EI a), to the positron- hydrogen s-wave annihilation parameter Z^b, 
calculated with the positron wavefunction in the static atomic potential (sol id curve), and with the Dyson orbital (dashed 
curve). Circles connected by dotted curve are the accurate results of Refs. [STf. 



The difficulty of calculating the vertex corrections to ZcS accurately is evident from these graphs, as all the diagrams 
in Fig. |5l contribute significantly. The higher-order diagrams containing the vertex function, are close to or greater 
than the lower-order diagrams. Note that all contributions have a similar dependence on the positron momentum. It 
is driven by the momentum dependence of the positron Dyson orbitals (external lines in the diagrams) , as discussed 
in Sec. ITlDl 

Figure ITTI shows that for p and d waves, the contribution of the diagram Fig. I^f ) grows rapidly and becomes largest 
towards the Ps formation threshold. This diagram describes annihilation inside the virtual Ps formation, which has a 
vigorous energy dependence close to threshold This may have a bearing on the kinks in Zgff visible at k ^ 0.6 a.u. 
for the p and d waves. Although they may be a numerical artifact, an indication of an inflection point is also present 
in the accurate p-wave results of Ref. |53| . 

The final results for i? = 15 a.u. are slightly higher than for i? = 30 a.u. Indeed, annihilation takes place near the 
atom, and a denser knot sequence for i? = 15 may provide a better description of small electron-positron separations. 
Our final values compare well with the accurate positron-hydrogen results obtained using a variational approach [53 • 
In that work the electron-positron distance was represented explicitly in the calculation, while we use a single-centre 
expansion. We believe that the remaining small discrepancy could be eliminated by 'pushing harder' the numerics in 
our approach. Thus, Fig. 1 121 shows that if we use n = 60 B-splines of order fc = 9 and include the first 23 basis states, 
the difference between our Zcs and the benchmark values is halved. 

The main conclusion of this section is that the numerical implementation of the many-body theory approach 
proposed in this paper works. For positron collisions with hydrogen, where this approach is exact, the calculations 
reproduce the best scattering phaseshifts, and yield good results in the more difficult annihilation problem. 



V. CORRELATION CONTRIBUTION TO THE ANNIHILATION VERTEX. 



The ability of a many-body theory to describe correlation corrections to the annihilation vertex may give some 
insight into the role such corrections play in positron annihilation with matter. Thus, in theoretical studies of positron 
annihilation in condensed matter, one usually starts by performing a calculation for the ground state electron and 
positron densities. The annihilation rate is then found via an equation of the form (see, e.g., [^), 



A = 7rr^c / pe(r)pp(r)7(pe,pp)dr. 



(37) 
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FIG. 11: Annihilation parameter Z^n for the s-, p-, and d-wave positron on hydrogen. Contributions of individual diagrams 
from Fig. |S]are: thick solid line, diagram (a); thick dotted line, (b); long dashed line, (c); thick dashed line, (d); dot-dashed 
line, (e); thin dotted line, (f) (all for f? = 15 a.u.). Thin solid line is the total for i? = 15 a.u., thin dashed line, total for 
7? = 30 a.u. Circles are the results of Ref. 153 . 



where /3p(r) and Pe(r) are the positron and electron densities and ^{pe, Pp) is the so-called enhancement factor. It is 
introduced to account for the Coulomb attraction in the annihilating electron-positron pair. It has long been known 
that the independent-particle approximation (7 = 1) underestimates the annihilation rates by several times j59j . and 
a number of semi-empirical and interpolation forms of ^{pe, Pp) have been suggested (see, e.g., Ref. |63| and references 
therein) . 

A comparison between Eq. (|37|l with 7 = 1 and the many-body theory expression l|24|) for the annihilation rate, 
shows that the former corresponds to the Oth-order term in Zcf^^ Eq. (|22|l [6l|. As we have seen in Sec. lIVI fFig. llOfl . 

does underestimate the annihilation rate in hydrogen by a factor of 5, even when the best positron wavefunction 
is used. On the other hand, the correlation correction to the annihilation rate [2nd term in Eq. H24|l ] does not have the 
form of Eq. H37|l . This correction depends on the positron wavefunction at two different points, r and r'. Therefore, 
any local expression like Eq. H37|l is necessarily an approximation. 

To illustrate this point, a contour plot in Fig. 1131 shows the radial part of the integrand of the nonlocal term. 



P,z(r)A«(r,r')Pei(^'), 



(38) 



for the s-wave (/ — 0) positron annihilation on hydrogen at fc = 0.06 a.u. The plot confirms that the correlation 
contribution to the annihilation vertex is localised near the atom. Its maximum at r = r' w 1.5 a.u. compares 
well with the radius of the hydrogen atom, and the ridge-like structure indicates that the "nonlocality" is limited to 
\r — r'l ~ 1 a.u. 
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FIG. 12: Positron annihilation on hydrogen, s wave. Many-body theory results: dashed line, 7? = 30 a.u. 
15 states used; dot-dashed line, R — 30 a.u., n — 60, k — 9, first 23 states used; circles, Ref. fs?]. 



n = 40, k = 6, first 



The overall size of the (nonlocal) correction to the annihilation vertex can be characterised by the average enhance- 
ment factor 7, 



7 



1 



/ V*(r)Ae(r,r')V'e(i-')t^rdr' 
IEn\^n{rW\MrWdr ' 



(39) 



We defined this factor in such a way that when used in place of 7 in Eq. H37|l . together with Pe(r) = ^'^d 
Pp(r) = |?/'£(r)P, it would reproduce correct values of the positron-atom annihilation rate. 

The quantity 7 can also be defined as the ratio of the total Zcs to the value obtained from the Oth-order diagram, 
Zj^g , Fig. El^a). Because of the weak energy dependence of (see Sec. Ill D|) . this ratio should not depend strongly 
on the energy of the incident positron. This also implies that 7 will be relatively insensitive to the wavefunction 
used to describe the incident positron. In particular, the use of either HF or Dyson wavcfunctions i/'e for the incident 
positron should yield close values of 7. 

As a test, this ratio was evaluated for the annihilation of the s, p and c?-wave positrons on hydrogen, see Fig. 1141 
As expected, the values of 7 depend weakly on the incident positron energy, except when the Ps formation threshold 
is approached. This confirms the earlier observation that the momentum dependence of various contributions to Z^g 
in Fig. 111! is approximately the same. The results obtained within the static approximation and with the Dyson 
orbitals are also very close to each other, even though the absolute values of Z^s obtained in the two approximations 
are very different (see, e.g.. Fig. IIUI) . One cannot help noticing that the enhancement due to annihilation vertex 
corrections increases with the angular momentum of the positron. For larger the penetration of the positron into 
the electron-rich regions of the atom is suppressed by the centrifugal barrier, and the effect of "pulling the electron 
out" (i.e., virtual Ps formation), described by the correlation corrections, becomes more important. 

Equation (|37|l is local, whereas equation H24|) contains a non-local term. In order to compare the enhancement of 
Zcs due to this term with the enhancement factors used in condensed matter calculations, we need to "localise" the 
contribution of A^{r, r'). Let us do this by introducing an effective 'electron correlation density' Pe{r) via the relation 



P,i{r)Ai'\r,r')P,i{r')drdr' = / p,{r)P^iir)dr, 



(40) 



for positron annihilation in the Ith partial wave. This allows us to define an effective enhancement factor ^e{r) through 



Pe{r)+Pe{r) = ■ye{r)pe{r), 



(41) 
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FIG. 13: Radial dependence of the integrand (1381 of the correlation correction for the annihilation of the s-wave positron with 
momentum k — 0.06 a.u. on hydrogen. The "ripples" is an artifact of the reconstruction of Ai'\r,r') from the matrix element 



which is equivalent to 

jeir) = l + Pe{r)/Pe{r), (42) 

where Pe{r) = Pf^ (r)/47rr^ , and Pis — 2re^^ for hydrogen. Equations H40|) and 142() guarantee that if we use jeif) in 
Eq. (|37l) . correct annihilation rates will be recovered. 

Equation H4UII does not define Peir) uniquely. We use two different methods to obtain it numerically. The first one 
states, 

/A?'(./)P.,(.-M/ 

Pei{r) 

which ensures that p\^\r) satisfies Eq. H4Q|I exactly. However, it has a disadvantage in that it may have unphysical 
poles at the zeroes of the positron wavefunction. A second method involves calculating 

AW(r + a/2,r-a/2)dt7, (44) 

-2r 

which follows from Eq. (|4()|l if we change variables r, r' to r ± a/2, and keep the lowest order term in the Taylor 
expansion of Peiir + <^/2)P^i{r — cr/2). This approximation may not be accurate for small r, where the positron 
wavefunction varies rapidly, but should be correct for larger r values, where the peaking of ■'(r, r') at r = r' (Fig. 
I13II means that Peiir) varies slowly on the scale of typical a. 

Figure ITKl shows both electron correlation densities calculated for the s-wave positron with momentum k = 0.5 a.u. 
Apart from the small range of distances near the origin, the values of pe from the two methods are close, although 
p^e^ shows some numerical "noise" related to the reconstruction of Ae''(r, r') from its matrix elements. A comparison 
with the hydrogen ground-state electron density shows that the latter drops much faster with the distance from the 
nucleus. In fact, Pe is much greater than p^ at those r where the positron density is large, in agreement with the 
correlation contribution to Z^s being five times Z^^^ . 
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FIG. 14: Ratio 7 of the total Z^a to for positron annihilation on hydrogen. Values obtained with positron wavefunction 
in the static approximation are shown by solid curves, and those with the Dyson orbitals, by dashed curves. 

0.4 



0.3 



3 



>^ 0.2 



CO 

CD 
Q 



0.1 




12 3 4 

Radius (a.u.) 

FIG. 15: Effective electron correlation densities p\^\r) (dashed curve) and p?'(r) (dot-dash curve) for the annihilation of the 
s-wave positron with k — 0.5 on hydrogen. Shown for comparison are the ground-state electron density in hydrogen Pe{r) (solid 
curve), and the positron density P^i(r) (dotted curve, arbitrarily scaled). 

The electron correlation densities obtained above allow us to calculate the corresponding enhancement factors 
and compare them with a parametrisation of 7(pe,Pp) derived by Arponen and Pajanne |63 | for a positron in a 
homogeneous electron gas, 



7AP = 1 + 1.23rs 



0.0742r^ + -rf , 
5 



(45) 
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where is a measure of the average distance between the electrons, 

1/3 



(46) 



This enhancement factor is plotted in Fig. 1161 together with the effective enhancement factors 7e(?') derived from 
Eq. (|42|) using p\^\ which is more stable numerically than Note that unlike 7ap, the factor 7e obtained from 
our many-body theory approach, depends on the energy and orbital angular momentum of the positron. The values 
shown in Fig. 1161 correspond to the momenta, k ~ 0.06 (s wave) and 0.5 a.u. (s, p and d waves). 
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FIG. 16: Enhancement factors for positron annihilation on hydrogen. Solid curve, 77IP, Eq. (I45II : dotted hne 7e for s-wave 
positrons at fc = 0.06; dashed Iine,7e for s-wave positrons at A; = 0.5; long-dashed line, 7e for p-wave positrons at fc = 0.5; 
dot-dashed line, 7e for d-wave positrons at A: = 0.5. 

A feature common to all enhancement factors in Fig. 1161 is their rapid rise with the distance from the nucleus. 
This increase is related to the drop of the electron density, a relation which is explicit in Eq. H45|) for 7ap . At small 
distances, where the electron density is large, 7ap compares well with 7e. However, at larger distances, where the 
electron density is low, 7ap is much greater than all of the 7e. Such a discrepancy could be expected, given that a 
homogeneous electron gas theory used to derive tap is more reliable in the high-density limit. The exaggeration of 
the enhancement by 7ap was reported in Ref. |6nl | where various forms of the enhancement factor were tested by 
comparison with accurate annihilation rates for a number of positron-atom bound systems. For Be and Mg (which 
have / > 6.8 eV), the values obtained using 7ap overestimated the accurate annihilation rates by factors of 5 and 
2, respectively. It must be noted that the product /0e7AP remains finite as pe —>■ 0. This means that one cannot in 
principle use it in Eq. H37I) for a continuous spectrum positron, since it would yield infinite values of the annihilation 
rate and Z^g. For the same reason a much stronger overestimate observed in Ref. !6(| with 7ap for Be is a direct 
consequence of the positron binding energy for Be being much smaller than for Mg. 

Figure [TBI confirms that 7e is energy-dependent. It is greater for k — 0.5 a.u. compared with k — 0.06 (s wave), the 
higher momentum being closer to the Ps formation threshold. Values of 7e derived for the p and d waves are similar 
to those from the s wave. Hence, the large differences between the average enhancement factors in Fig. El are due to 
the effect of the centrifugal barrier on the positron wavefunction. The distances which effectively contribute to Zcs 
are greater for the positron in higher partial waves. 
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VI. SUMMARY AND OUTLOOK 



In this paper we have formulated a many-body theory approach which accounts for the main correlation effects 
in positron-atom interactions. These are (i) polarisation of the target by the positron, (ii) virtual positronium 
formation, and (iii) strong enhancement of the electron-positron contact density due to their Coulomb interaction. 
The key development for an accurate description of (ii) and (iii) is the summation of the ladder diagram series and 
calculation of the electron-positron vertex function. B-spline basis sets and extrapolation over the orbital angular 
momenta are used to achieve convergence of the sums over the electron and positron intermediate states. The method 
can be applied to a range of problems such as positron scattering, annihilation and formation of bound states. 

Although our main interest is in exploring many-electron targets, the method has been first tested for hydrogen, 
where accurate benchmark data exist for the scattering phaseshifts and annihilation rates. In the case of hydrogen 
the present formalism is exact. Numerically, excellent agreement with accurate variational calculations for the phase 
shifts has been obtained, together with a good agreement for the annihilation parameter ZcS- The calculation of 
the most difficult part of the correlation potential, which contains the vertex function, for many-electron atoms is 
only marginally more difficult than for hydrogen. Therefore, we expect that application of our many-body theory to 
the problems of positron scattering and annihilation on noble-gas atoms and binding to halogen ions [s^, 0| should 
yield accurate results. In particular, we would like to re-examine and improve the accuracy of the many-body theory 
predictions [s^ of positron binding energies to the ns"^ atoms such as Mg, Cd and Zn. 

The advantage of many-body theory methods is their physical transparency. It allows one to distinguish between 
different physical mechanisms and compare their relative importance. Thus, we saw that virtual Ps formation in 
positron-hydrogen scattering is just as important as the target polarisation. Correlation corrections to the annihila- 
tion vertex, which are physically related to the virtual Ps formation, are even more important. They enhance the 
annihilation rate in hydrogen by a factor of 5 or more, depending on the positron partial wave. Such vertex corrections 
depend weakly on the positron energy, and the enhancement they produce is practically the same for various positron 
wavefunctions. Therefore, one could use the average enhancement factors derived for isolated atoms to obtain reli- 
able annihilation rates for atoms placed in different environments, provided that a sufhciently accurate single-particle 
positron wavefunction is available. Of course, different atomic subshells will be characterised by different enhancement 
factors. However, they can all be determined in an many-body calculations of the type described in this paper, and 
serve as an input for the calculations of positron annihilation in molecules or condensed matter. 

In this work we have analysed the spatial dependence of the nonlocal correlation corrections to the annihilation 
vertex. We have also derived the equivalent local enhancement factor and compared it with an expression used 
in condensed matter calculations. Similar comparisons for larger many-electron targets may test various forms of 
enhancement factors in a much greater range of electron densities. 

The rapid development of computers over the past few decades seems to have favoured theoretical methods other 
than many-body theory. Such methods often rely more on the computer power and numerical techniques than on the 
physical insight. They often appear to be "more exact" than the sophisticated but explicitly approximate many-body 
theory approaches, and promise improved results due to shear growth of computer power. Their drawback is that they 
do not always increase one's physical understanding of the problem. We believe that a further theoretical development 
of many-body methods combined with a judicial use of computers is a healthy alternative. 
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For a positron interacting with a spherically-symmetric target, both the self energy S£;(r,r') and the correlation 
correction to the annihilation vertex A£(r,r') can be expanded in partial waves, e.g., 
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APPENDIX A: TECHNICAL DETAILS 




(Al) 



A=0 f_i=-\ 
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Defining the volume element dr = r'^drdfl, and the positron wavefunction with orbital angular momentum /, (psi^) 
r^^ P^i{r)Yim{^), we obtain the matrix element (e|S|e') as, 



I ^:(r)I]B(r,r')^e'(r')rfr 



= P,i{r)Y!^^{r,r')P,,i{r')drdr'. (A2) 



{e\^E\e') 



The angular reduction of the various diagrams in (£|I]£;|e'} and {e\/S.E\£) is simplified by the use of graphical techniques 
for performing angular momentum algebra |65l |. The final expressions for the diagrams in terms of reduced matrix 
elements, 3j and 6j symbols are given below. 

The reduced Coulomb matrix element is defined as, 

ooo/looo 



(3,4||FH|2,1) = VMMMh] 

X j P,^i^{n)P,,i,{r2)^Pe,i,ir2)PeMri)dridr2, (A3) 

where [^i] = 2/i + 1, etc. Note that instead of including the minus sign in the reduced Coulomb matrix element when 
it involves the positron, we account for it in the overall sign factor for the diagram (see below). The reduced Coulomb 
matrix element for an electron-positron pair coupled into a total angular momentum J is given by 

(3,4||F(^)||2,1) = ^(-l)-^+'(3,4||;^||2,l)|"[|^|^|. (A4) 

This expression is similar to the 'exchange' matrix element that one meets in all-electron problems |22l |. 

The sum of the ladder diagram series (the vertex function) is calculated via the matrix equation © , which must 
be solved for all possible total angular momenta J of the electron-positron pair: 

/ iiT-(-')ii \ / \\^dJ)\\ \ >^ {^2, ^J.2\\v^•'^\^^,^){v, ^i\\T^E\\^^l,^l) 
{v2,^l2\\T'■E'Ul,vl) = -{y2,^i2\\v'~''>\\^Ml,vl) _ . (A5) 

In the Appendix, the state labels v, vi, etc., refer to the positron orbitals e^l^, £viKi, etc. Similarly, ^, etc., 
label excited-state electron orbitals. Electron orbitals occupied in the target ground state (holes) are labeled by Latin 
indices (n). When B-spline basis states are used, Eq. (|A5|) is a finite-dimension matrix equation solved by matrix 
inversion. 

The self-energy diagrams can be expressed in terms of these matrix elements. For closed-shell atoms, each loop in 
the diagram gives a spin factor of 2. This factor should be omitted for hydrogen which has only one electron in the 
Is orbital. The sign factor for each diagram is (— l)""^'''^'^, where a is the number of hole lines, h is the number of 
electron-hole loops and c is the number of positron-electron Coulomb interactions. In the expressions below we label 
the incident positron angular momentum by Ip. This is also the orbital angular momentum of the external lines e, e' 
of the diagrams. 

The 2nd-order self-energy diagram. Fig. ^a), is given by 

The virtual-Ps contribution to E^;, Fig. |21 is obtained after finding the vertex function as follows: 

2 ^ ^ [J](£^ n\\V(J^ t^2)(^2, M2||ri;^L \\fii,i^i){i^i,fii\\V^'^ s) ^^^^ 

Matrix elements of the annihilation (5-function are defined similarly to the Coulomb ones, using the expansion of 
the (5(ri — r2) in terms of spherical harmonics The reduced matrix element then is 



X / P,^i^ir)P,,uir)P,,i,{r)P,,i,{r)r-'dr, (A8) 
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and the matrix element for an electron-positron pair coupled into a total angular momentum J is given by, 

(3,4||5('^)||2,1)=5](-1)''+'(3,4||5,||2,1)|'[J^J^|. (A9) 

The various diagrams contributing to Z^s can then be expressed in terms of these matrix elements. As with the 
self energy diagrams, the factor of 2 which results from closed loops, must be removed for hydrogen. Note that in the 
diagrams below, we also include factors of 2 to account for the mirror images of those diagrams that are not symmetric. 
To produce correct values of Z^s the expressions given below must be multiplied by the positron normalisation factor 

The Oth-order diagram, Fig. 01a), is a sum of simple radial integrals over all hole orbitals n, 

2ET^/^elMnl.„(-K^rf- (AlO) 
The Ist-order contribution, Fig. Efb) and (c), is given by 

[l][lp]{E + e„ -e„ -e^) 



Expressions for the remaining four contributions, Fig. [Itc)-(f), are: 



2 ^ J2 W(£,^r''^^llM2,^2)(^2,M2p(^'Hm,'^i)('^i,Mi||^^-^^ll»,£) ^ (^^2) 



[J] (e, ||M2, ^^2)(i'2, M2|l4^|,Jln, e) 

Vi,iJii,n J 



. S-^ ST^ ['^J\t, ' 11/^2, '^2/\^2,ft2||^g+e„ir-,^/ (AT-i) 



[J](£,n||V^('^)||M3,^^3)(i^3,M3||'5('^)||A*2,^^2)(i'2,M2|i4+eJI".e) 



(AM) 



2 [J] "P^Ie^ IIM3, t^3)('^3, Mall'^^"^^ IIM2, ^2) (i^2, M2Pyjg,. H», s) ^^^^^ 

^. „ J [^p\{E -\- £n — Ev3 ~ £/i3)(-E' + £n ~ £1^2 ~ £^2) 

where we have introduced 

{v2,M\A'eUM.^) ^ 2^ _ _ (A16) 

This amplitude not only helps to make the expressions above more compact. It also provides a practical way of 
computing these diagrams. Thus, without a preliminary evaluation of {v2t ^J'■2\\A^E+g ||'^i£)i diagram (f) in Fig. [S] 
contains eight-fold summation over hundreds of electron and positron basis states. 



[1] W. J. Cody, J. Lawson, H. S. W. Massey, and K. Smith, Proc. Roy. Soc. A 278, 479 (1964). 

[2] M. Ya. Amusia, N. A. Cherepkov, L. V. Chernysheva and S. G. Shapiro, J. Phys. B 9, L531 (1976). 

[3] V. A. Dzuba, V. V. Flambaum, W. A. King, B. N. Miller, and O. P. Sushkov, Phys. Scr. T46, 248 (1993). 

[4] G. F. Gribakin and W. A. King, J. Phys. B 27, 2639 (1994). 

[5] New Directions in Antimatter Chemistry and Physics, edited by C. M. Surko and F. A. Gianturco (Kluwer Academic 

Publishers, Dordrecht, 2001). 

[6] M. J. Puska and R. M. Nieminen, Rev. Mod. Phys. 66, 841 (1994). 

[7] A. David, G. Kogel, P. Sperr and W. Triftshauser, Phys. Rev. Lett. 87, 067402 (2001). 



24 



[8] G. Jerusalem, R. Hustinx, Y. Beguin and G. Fillet, Eur. J. Cancer 39, 1525 (2003). 

[9] C. Schwartz, Phys. Rev. 124, 1468 (1961). 
[10] A. K. Bhatia, A. Temkin, R. J. Drachman, and H. Eiserike, Phys. Rev. A 3, 1328 (1971); A. K. Bhatia, A. Temkin, and 

H. Eiserike Phys. Rev. A 9, 219 (1974). 
[11] J. W. Humberston, J. Phys. B, 6, L305 (1973); R. I. Campeanu and J. W. Humberston, ibid. 10, L153 (1977). 
[12] J. W. Humberston, Adv. At. Mol. Phys. 15, 101 (1979). 

[13] J. W. Humberston, P. Van Reeth, M. S. T. Watts and W. E. Meyerhof, J. Phys. B 30, 2477 (1997). 
[14] P. Van Reeth and J. W. Humberston, J. Phys. B 32, 3651 (1999). 

[15] M. T. McAlinden, A. A. Kernoghan, and H. R. J. Walters, Hyperfine Interactions 89, 161 (1994); J. Phys. B 29, 555 

(1996); 29, 3971 (1996); 30, 1543 (1997). 
[16] G. G. Ryzhikh, J. Mitroy, and K. Varga, J. Phys. B 31, 3965 (1998). 

[17] V. A. Dzuba, V. V. Flambaum, G. F. Gribakin and C. Harabati, Phys. Rev. A 60, 3641 (1999); V. A. Dzuba, V. V. Flam- 

baum and C. Harabati, Phys. Rev. A 62, 042504 (2000). 
[18] I. A. Ivanov, J. Mitroy, and K. Varga, Phys. Rev. Lett. 87, 063201 (2001); J. Mitroy and I. A. Ivanov, Phys. Rev. A 65, 

012509 (2002). 

[19] M. W. J. Bromley and J. Mitroy, Phys. Rev. A 65, 062505 (2002); 65, 062506 (2002). 
[20] M. W. J. Bromley and J. Mitroy, Phys. Rev. A 66, 062504 (2002); 67, 062709 (2003). 

[21] A. L. Fetter and J. D. Walecka, Quantum Theory of Many-Particle Systems (McGraw-Hill, New York, 1971). 
[22] M. Ya. Amusia and N. A. Cherepkov, Case Studies in Atomic Physics 5, 47 (1975). 
[23] H. P. Kelly, Phys. Rev. 131, 684 (1963). 

[24] M. Ya. Amusia, N. A. Cherepkov, L. V. Chernysheva, S. G. Shapiro, and A. Tanchich, Zh. Eksp. Teor. Fiz. 68, 2023 (1975) 

[Sov. Phys. JETP 41 1012 (1975)]. 
[25] M. Ya. Amusia, N. A. Cherepkov, L. V. Chernysheva, D. M. Davidovic and V. Radojevic, Phys. Rev. A 25, 219 (1982). 

According to this paper the coefficient in Eq. 11511 is — vr, because they express the self-energy operator E_b in Rydberg. 
[26] W. R. Johnson and C. Guet, Phys. Rev. A 49, 1041 (1994); 64, 019901 (2001). 

[27] L. V. Chernysheva, G. F. Gribakin, V. K. Ivanov and M. Yu. Kuchiev, J. Phys. B 21, L419 (1988); G. F. Gribakin, 

B. V. Gultsev, V. K. Ivanov and M. Yu. Kuchiev, J. Phys. B 23, 4505 (1990). 
[28] W. R. Johnson, J. Sapirstein, and S. A. Blundell, J. Phys. B 22, 2341 (1989). 

[29] V. A. Dzuba, V. V. Flambaum, G. F. Gribakin and O. P. Sushkov, Phys. Rev. A 44, 2823 (1991). 
[30] V. A. Dzuba and G. F. Gribakin, Phys. Rev. A 49, 2483 (1994); 50, 3551 (1994). 

[31] V. A. Dzuba, V. V. Flambaum, P. G. Silvestrov and O. P. Sushkov, J. Phys. B 20, 1399 (1987); Phys. Lett. A 131, 461 

(1988); V. A. Dzuba, V. V. Flambaum and O. P. Sushkov, Phys. Lett. A 140, 493 (1989). 
[32] S. A. Blundell, W. R. Johnson and J. Sapirstein, Phys. Rev. A 38, 4961 (1988). 
[33] V. A. Dzuba, V. V. Flambaum, G. F. Gribakin and W. A. King, Phys. Rev. A 52, 4541 (1995). 
[34] V. A. Dzuba, V. V. Flambaum, G. F. Gribakin and W. A. King, J. Phys. B 29, 3151 (1996). 
[35] G. F. Gribakin and W. A. King, Can. J. Phys. 74, 449 (1996). 
[36] G. G. Ryzhikh and J. Mitroy, Phys. Rev. Lett. 79, 4124 (1997). 

[37] G. G. Ryzhikh and J. Mitroy, J. Phys. B 31, L401 (1998); J. Mitroy and G. G. Ryzhikh, ibid. 32, 1375 (1999). 
[38] J. Mitroy, M. W. J. Bromley and G. G. Ryzhikh, J. Phys. B 35, R81 (2002). 
[39] J. Mitroy, M. W. J. Bromley and G. G. Ryzhikh, in Ref. pp. 199-221. 

[40] A. B. Migdal, Theory of Finite Fermi-Systems and Applications to Atomic Nuclei (Interscience, New York, 1967). 
[41] J. S. Bell and E. J. Squires, Phys. Rev. Lett. 3, 96 (1959). 

[42] L. D. Landau and E. M. Lifshitz, Quantum Mechanics, 3rd ed. (Pergamon, Oxford, 1977). 
[43] For electrons this also means the absence of the contributions of the target exchange potential. 

[44] RPA provides reliable values of atomic dipole polarisabilities, see, e.g., M. Ya. Amus'ya, N. A. Cherepkov, and S. G. Shapiro, 
Zh. Eksp. Teor. Fiz. 63, 889 (1972) [Sov. Phys. JETP 36, 468 (1973)]; D. Kolb, W. R. Johnson, and P. Shorer, Phys. Rev. 
A 26, 19 (1982); W. R. Johnson, D. Kolb, and K.-N. Huang, At. Data and Nucl. Data Tables 28, 333 (1983). 

[45] G. F. Gribakin, in Photonic, Electronic and Atomic Collisions (XXII ICPEAC), Proceedings, Eds. J. Burgdorfer, J. Cohen, 
S. Datz, and C. Vane, (Rinton Press, Princeton, 2002), pp. 353-64. 

[46] P. A. Eraser, Adv. At. Mol. Phys. 4, 63 (1968). 

[47] This cross section is given by CT27 = nrQlc/v), where v is the electron-positron relative velocity, see, e.g., V. B. Berestetskii, 
E. M. Lifshitz, and L. P. Pitaevskii, Quantum electrodynamics (Pergamon, Oxford, 1982). Annihilation into two gamma 
quanta is allowed only for the zero total spin of the pair, S = 0. For 5 = 1 they annihilate into three photons. The 
corresponding spin-averaged cross section is almost 400 times smaller: a^j = [4(7r^ — 9)/3]aro(c/u), where a — e^ /he ~ 
1/137. 

[48] G. F. Gribakin and J. Ludlow, Phys. Rev. Lett. 88, 163202 (2002); J. Ludlow and G. F. Gribakin, Phys. Rev. A 66, 064704 
(2002). 

[49] J. Sapirstein and W. R. Johnson, J. Phys. B 29, 5213 (1996). 

[50] C. de Boor, A Practical Guide to Splines (Springer, New York, 1978). 

[51] M. Ya. Amusia, L. V. Chernysheva, Computation of Atomic Processes. A Handbook for the ATOM Programs (lOP Pub- 
lishing, Bristol, 1997). 

[52] I. Bray and A. T. Stelbovics, Comp. Phys. Comm. 85, 1 (1995); I. Bray and D. V. Fursa, Phys. Rev. Lett. 76, 2674 (1996). 
[53] I. Bray and A. T. Stelbovics, Phys. Rev. A 48, 4787 (1993). 
[54] J Mitroy and G Ryzhikh, J. Phys. B 32, 2831 (1999). 



25 



[55] G. F. Gribakiri and J. Ludlow, J. Phys. B 35, 339 (2002). 

[56] C. Schwartz, Phys. Rev. 126, 1015 (1962); Methods in Computational Physics, Vol. 2 (Academic, New York, 1962), pp. 
241-266. 

[57] P. Van Reeth and J. W. Humberston, J. Phys. B 31, L231 (1998). 

[58] G. F. Gribakin, Phys. Rev. A 61, 022720 (2000). 

[59] R. A. Ferrell, Rev. Mod. Phys. 28, 308 (1956). 

[60] J. Mitroy and B. Barbiellini, Phys. Rev. B 65, 235103 (2002). 

[61] This is true to the extent that we neglect the small difference between the true and HF electron densities, and that between 

|i/>^(r)|2 and pp(r). 

[62] J. Arponen and E. Pajanne, Ann. Phys. 121 343 (1979). See also B. Barbiellini, M. J. Puska, T. Torsti, and R. M. 

Nieminen, Phys. Rev. B 51, 7341 (1986). 
[63] J. Ludlow and G. F. Gribakin, to be submitted to Phys. Rev. A (2004). 
[64] J. Ludlow and G. F. Gribakin, to be submitted to J. Phys. B (2004). 

[65] D. A. Varshalovich, A. N. Moskalev and V. K. Khersonskii, Quantum Theory of Angular Momentum (World Scientific, 
Singapore, 1988). 



